Microscopic chaos from Brownian motion in a ID anharmonic oscillator chain 
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The problem of relating microscopic chaos to macroscopic behavior in a many-degrees-of-freedom 
system is numerically investigated by analyzing statistical properties associated to the position and 
momentum of a heavy impurity embedded in a chain of nearest-neighbor anharmonic Fermi-Pasta- 
Ulam oscillators. For this model we have found that the behavior of the relaxation time of the 
momentum autocorrelation function of the impurity is different depending on the dynamical regime 
(either regular or chaotic) of the lattice. 
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I. INTRODUCTION 

A physically relevant and still not satisfactorily solved 
problem is the question of how the underlying chaotic 
microscopic dynamics of many-particle systems is related 
to the observed macroscopic behavior. Assuming that 
the main features observed at the macroscopic level can 
be accounted for by the microscopic dynamics, one would 
expect that some characteristics of the latter, specifically 
its chaoticity, should be detected, in principle, at the 
macroscopic level of description. 

Recent empirical evidence suggesting microscopic 
chaos on a molecular scale has been presented in an ex- 
periment on the position of a Brownian particle (BP) 
in a fluid The measurements were made at regu- 
lar time intervals and the experimental time series data 
was then interpreted using standard techniques of chaotic 
time series analysis, suggesting a positive lower bound on 
the Kolmogorov-Sinai entropy, hence, microscopic chaos. 
However, a similar bound has been obtained with com- 
puter experiments on the nonchaotic Eherenfest wind- 
tree model where a single particle diffuses in a plane due 
to collisions with randomly placed, fixed, oriented square 
scatterers , rendering doubts about the conclusion that 
microscopic chaos has been experimentaly detected. Fur- 
ther comparisons with the chaotic Lorentz model, which 
has circular scatterers and has a diffusive behavior that is 
undistinguishable from the one exhibited by the Eheren- 
fest model, have confirmed that the standard methods of 
chaotic time series analysis are ill suited to the problem 
of distinguishing between chaotic and nonchaotic micro- 
scopic dynamics, although some alternatives have been 
proposed in order to overcome this situation . 

The results in Ref. ^ were attributed to the physical 
issue of time scales - the time interval Ai « 1/60 s be- 
tween measurements was vastly greater than the typical 
collision times tc ~ 10^^^ s of the BP with the solvent 
particles in the fluid, which make the diffusive behavior in 



the experimental data and in the computer experiments 
virtually identical ||^. However, there is a time scale in 
Brownian motion, not mentioned in the previous works, 
that is characterized by the relaxation time r « 10~* s 
of the momentum autocorrelation function (MACF) of 
the BP. Since At ^ t > tc, this relaxation time could, 
in principle, be used to probe the microscopic dynamics 
of the system, but so far this possibility has not been 
explored experimentally. Moreover, the simplicity of the 
models used in |^,^ precludes a study in the momentum 
space from the very beginning, and therefore, a different 
kind of model is needed to address the problem at hand. 

We propose to reconsider a simple Hamiltonian model 
that has been studied extensively in the past, namely, a 
one-dimensional chain of harmonic oscillators of unitary 
mass coupled to a heavy oscillator (impurity). Many- 
body effects in Brownian motion can be analytically in- 
vestigated with this model, and it was proved that the 
impurity of mass Af 3> 1 satisfies a Langevin equation 
when the thermodynamic limit — > oo is taken and ran- 
dom initial conditions are given to all the oscillators of 
the system [|]. However, in its original formulation this 
model is useless to address the proposed problem because 
of the lack of chaoticity in the lattice, which nevertheless 
behaves as a heat bath. 

In this paper we extend the above-mentioned harmonic 
model with an anharmonic potential, which makes the 
resulting Hamiltonian very similar to the Fermi-Pasta- 
Ulam (FPU) model With this modification we ex- 
plore, for this particular model, the possible effect of mi- 
croscopic chaos in the statistical behavior of the heavy 
impurity. We note that the problem of an harmonic os- 
cillator coupled linearly to a FPU chain has already been 
studied Q . Diatomic and disordered mass versions of the 
FPU model have been applied to the problem of energy 
equipartition ||^ and heat conduction but not with a 
single heavy impurity in the context of Brownian motion, 
as is our case. 
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II. THE MODEL 

The Hamiltonian of the model we are considering can 
be written, in terms of dimensionless variables, as 
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where = 1 if i ^ and mo — M; periodic bound- 
ary conditions are assumed (a;(jv/2)+i = The 
model describes a system of one-dimensional N coupled 
nonlinear oscillators of unit mass with nearest-neighbor 
interactions and a central oscillator (impurity) of mass 
M, with displacement xq = X and momentum pq = P. 
The quadratic part of Eq. was the one studied in 
connection with Brownian motion in Ref. Q. The value 
/3 = 0.1 was used in all the numerical experiments here- 
after reported. We call this model the modified FPU 
(MFPU) model. 



III. DYNAMICS IN PHASE SPACE 

It is known that, for high values of the total energy 
per degree of freedom e, the homogeneous (uniform mass) 
FPU model [obtained from Eq. (|^) by taking M = 1] is 
chaotic whereas, for small e values, the model behaves as 
a chain of harmonic oscillators |^ . It is then important 
to corroborate if these dynamical regimes are in some 
way affected by the presence of the impurity at lattice 
site i = Q. To this end we employ the following distance 
in phase space 

where the sum runs over all the iV -|- 1 oscillators of 
the system, the superscripts 1 and 2 refer to two states 
that at time t = differ by an infinitesimal quantity 
d(0) = 10^^ or less. In our simulations we prepared a 
system with N = 300, 000 unit mass oscillators in their 
equilibrium positions. Then we distributed the momenta 
of the oscillators according to a Maxwell-Boltzmann dis- 
tribution at a temperature T consistent with the chosen 
value of e and let the system evolve in time by solving 
the equations of motion with an improved leap-frog algo- 
rithm. All computations were carried out in double pre- 
cision. At each temperature an average ((d(t))) over ten 
different realizations was taken to avoid fluctuations that 
come from a particular choice of the initial conditions. 
This calculation was performed for different values of e 
and M = 1, 40, 60, 80, and 100. The results for M = 1 
(shown for comparison) and 60 are exhibited in Fig. ^. 
As can be seen, an exponential divergence corresponds 



to values e > 1 whereas for e < 1 the distance is found to 
be bounded within the same order of magnitude. This 
happens independently of the value of M for all the cases 
studied. Thus we can conclude that the qualitative dy- 
namical behavior of the system is only weakly affected by 
the presence of the heavy impurity. Hence, for small val- 
ues of e, which correspond to the value range 0.01 < e < 1 
in Fig. |l|, we find regular (i.e., almost periodic) behavior 
in phase space whereas for large e values, which we have 
studied in the range 1 < e < 10, the dynamics in phase 



space is chaotic |10 
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FIG. 1. Average distance in phase space between pairs 
of trajectories corresponding to initial conditions such that 
d{0) < IQ-^ and e = 0.01 (circle), 0.1 (square), 1 (plus), 
2 (cross), and 10 (asterisk) for M = 1 and M = 60 with 
N = 300, 000. Time is measured in natural units. 



IV. STATISTICAL BEHAVIOR OF THE HEAVY 
IMPURITY 

Thermal equilibrium between the impurity and the 
FPU chain was attained within the time scale t = 5x 10^. 
Time averages {■ ■ ■)t necessary to obtain the relevant 
correlation functions were computed afterwards over the 
characteristic time scale t = 2 x 10^. The first quantity 
to be studied was the mean-square displacement. An ex- 
ample is shown in Fig. ^; the temperature (i.e., mean 
kinetic energy per degree of freedom) T = 11.61 was 
computed, which is the same value reported in Ref. |l^] 
for the homogeneous FPU model and e = 10. For short 
times {t < 100) the ballistic behavior of the impurity is 
evident, whereas, for long times, the diffusive effect of the 
bath of light oscillators has been established. The noise 
in the curves reported is negligible for the time interval 
100 < t < 300, so we can be confident of the values of 
the "diffusion coefficients" Dj^^ obtained from the slope 
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FIG. 2. Mean-square displacement of the impurity for 
e = 10. The solid line, with slope 2D = 7.95, has been added 
for reference. 

An important issue at this point is to check if the devi- 
ations of the individual D values from their arithmetic 
mean D are statisticaly significant for any given e; if so, 
then there is a dependence of the diffusion coefficient on 
the mass of the heavy oscillator. For a system of coupled 
harmonic oscillators it is well known that the diffusion 
coefficient is indeed independent of the mass of the impu- 
rity. However, a different situation arises in another type 
of models. In Ref. |12| it was shown that for a system of 
noninteracting mass points moving on a line and collid- 
ing with a heavy particle there is a mass dependence that 
increases with the mass of the heavy test particle. For 
the MFPU model our results show that the deviations 
Dj^j — D are negligible for small e values and more pro- 
nounced for large e values. Nevertheless, in both regimes 
(e < 1 and e > 1) the deviations never exceeded a few 
percent in all the cases considered, with no systematic 
drift. Thus we can consider that the diffusion coefficient 
is given hy D ^ D. For e — 0.01 we have obtained the 
relationship 2D ~ 0.0104 k, T, which corresponds to the 
exact result for the harmonic chain Q. In Fig. | we 
plot the value of the diffusion coefficient D vs the energy 
density e. By a least-squares fit we obtain a power-law 
scaling with respect to e of the form D = -Doe™, where 
TO = 0.964 ± 0.008 and Dq = 0.466 ± 0.008. From the 
validity of this relationship throughout the entire e value 
range studied (0.01 < e < 10), from regular to chaotic, 
we conclude that, as far as our results allow us to infer, 
no evidence of the microscopic dynamics can be detected 
in the behavior of the diffusion coefficient, in agreement 
with the known results of the Eherenfest and Lorentz 
models ijl. 
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FIG. 3. Diffusion coefficient vs energy density. The 
straight line is the power-law fit explained in the text. 

The MACF po{t) = {P{t)P{Q))tl {P^{Q))t of the im- 
purity can also be straightforwardly computed. Figure ^ 
shows the time dependence of po{t) in linear- logarithmic 
scale for all values of M considered and e = 10; other e 
values yield similar results. It is evident that, for the rel- 
evant time scale (0 < t < 50), the MACF exhibits an al- 
most exponential decay, which allows us to obtain the re- 
laxation time T from the slope of each curve. The staight 
lines plotted correspond to the exponential fit exp(— i/r). 
The observed accuracy of this fit is a clear signature that 
the impurity indeed performs Brownian motion. The be- 
havior of the MACF is smoother as the mass M increases, 
so this fact can be considered as another evidence that 
the more physically relevant results are obtained with 
larger values of M . The magnitude of po{t) is negligible 
for t > 50 in all cases, and consequently its contribution 
was not considered in the computation of t. 
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FIG. 4. Logarithm of the MACF of the impurity for e = 10. 
Symbols have the same meaning as in Fig. ^. The straight 
lines for each value of M correspond to the fit exp(— t/r). 

A very interesting result can be obtained if we plot the 
relaxation time r vs the energy density e for the entire e 
value range, from regular to chaotic, and all the consid- 
ered M values, as shown in Fig. || in log-log scale. We 
observe that the data points are separated into two dif- 
ferent and well-defined regions, depending on the e value. 
In all cases the dependence of t on the energy density is 
weak when e < 1. On the contrary, when e > 1^ t de- 
creases steeply with increasing e. In each of these regimes 
T has a power-law scaling with respect to e as 
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T„„e"M, 0.01 <e<l 
e"Af , 1 < e < 10 . 



(3) 



For all the data points corresponding to a fixed M value, 
Tg 7^ „ and a^j ^ . This feature contrasts with 
the behavior of D as seen in Fig. ^, where a single scaling 
law is valid for the entire e value range. In the e < 1 
regime the accuracy of the power-law scaling of r^^ can be 
tested by noticing that, for e — 0.01, t^, must coincide 
with the exact value M/2 of the relaxation time of the 
MACF for the harmonic chain [Q. This is indeed the 
case with an accuracy better than 1%. 
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FIG. 5. Relaxation time vs energy density. The dashed 
(dotted) lines correspond to the least-squares power-law fit (^) 
for the e < 1 (e > 1) regime and all M values. 

We notice that the slopes of each of the fits in Fig. |^ 
are approximately the same for the e < 1 regime. Fur- 
thermore the same thing happens for e > 1, though with 
other slope value. These facts imply that the scaling ex- 
ponents in Eq. (jsj) for each regime are approximately the 
same, independent of the particular M value to which 
they correspond. It is then reasonable to expect the ex- 
istence of a mass-independent scaling exponent a* (a*) 
for the e < 1 (e > 1) regime. In order to compute a* (a*) 
we rescale the data of Fig. ^ on e < 1 (e > 1) with their 
corresponding r„ j,,, ) value to obtain the new data 

e* = In e/ InT^ (e* = Ine/ Inr^ ) and r* = Inr/ Inr^ 
(r* = Inr/lnTg ^j). The result is reported in Fig. ^. 
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FIG. 6. Same data as in Fig. ^, but in rescaled units r* 
and e*. Dashed lines are the least-squares fit in each region. 



As can be seen, we have found a common scaling for 
all the data points in each regime, e < 1 and e > 1. The 
dashed lines correspond to the best fit of the data in each 
region, which allow us to calculate the mass-independent 
scaling exponents a* and 5* . The result is 
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^g*^-0.017±0.004 
^g*^-0.136±0.007 



-1.57541 < e* < e* 
e! < e* < 0.80121. 



(4) 



e* « is the critical value at which the change of scaling 
occurs and is indicated by the dotted vertical line. The 
change in behavior when going from low to high e value 
is quite pronounced since the scaling exponents differ by 
one order of magnitude. 



V. DISCUSSION 

After the presentation of our results we would want to 
make some remarks. First of all, it could be argued that 
the diffusion coefficient can be obtained from the MACF 
by employing the Green-Kubo relation [l^ ] that can be 
written, in terms of the MACF, as O 
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(5) 



This identity implies that D and r are not independent 
variables. More precisely, if we substitute the exponential 
fit pa{t) — exp(— t/r) in Eq. (^) we obtain the relation 
Dj^^ = {T/M)t, and therefore, the behavior of r and 
Dj^j as e increases should be similar. However, we have 
T/M < 1 for aU M values considered, e.g., T/M = 2.5 x 
10"'' for e = 0.01 and 0.29045 for e = 10 with M = 40. 
Thus the behavior of e vs r depicted in Fig. |5|is lost when 
T is multiplied by the factor T/M to obtain 1?^,, . In fact, 
if we rescale each of the data points displayed in Fig. |^ 
with the corresponding factor T/M, the result are data 
sets {Dj^j}, one for each M value, spanning the entire e 
value range. Furthermore, if we take, for a certain e, the 
arithmetic mean of all the Dj^^ values, we recover the data 
points displayed in Fig. |^. Thus there is no contradiction 
between the data in Figs. || and |^ and the Green-Kubo 
relation (H). 

To address the problem of the origin of the change in 
the scaling of r as depicted in Fig. || we consider the 
dynamical behavior of the system as described in Fig. 
|l|. We notice that, for the particular value of M = 60 
(the following observations are valid for other M values 
as well), when e < 1 the system is in a regular region of 
phase space, whereas, for e > 1, the dynamics is chaotic. 
This e value range is approximately the same as that ob- 
tained from Fig. |^ for the change in behavior of r as a 
function of e. From this observation we conclude that the 
change of scaling in r may be attributed to the qualita- 
tive change, from regular to chaotic, in the microscopic 
dynamics of the whole system. 
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It is important to point out that similar numerical 
results were obtained from the behavior of the relax- 
ation time T„ of the spectral entropy, a function used to 
study the relaxation to equilibrium of a many-degrees- 
of-freedom Hamiltonian system starting with far-from- 
equilibrium initial conditions. For high values of e (total 
energy per degree of freedom) ~ const, but sud- 
denly increases below a certain threshold, ^ 1 for the 
homogeneous FPU model 1^. This result bears some re- 
semblance to the change of scaling in r for our model, al- 
though the details of behavior as e decreases are different 
for both relaxation times r and r^. Another difference 
worth stressing is that the relaxation time for the ho- 
mogeneous FPU model is associated to the equilibration 
process, but the relaxation time r reported in Fig. || is 
the time that characterizes the decay of the MACF of our 
MFPU model, which is computed in the thermodynamic 
equilibrium state . Moreover, it should be noted that 
numerical analysis of the dynamics of the homogeneous 
FPU Q and other one-dimensional models 1 16| has shown 
that the change of behavior in as e increases is related 
to a transition, common to these models, from weak (i.e., 
almost periodic) to strongly chaotic regimes. It remains 
an open problem whether the behavior of r as reported 
in Fig. ^ can be obtained or not when the impurity is 
coupled to some other system, such as those studied in 
Ref. ijl^. The sudden change in the scaling of r clearly 
observed in Fig. |^ could very well be related to some par- 
ticular characteristic of the aforementioned transition in 
the case of the FPU chain. We will address this problem 
in a future communication. 



VI. CONCLUSIONS 

We have found, for the particular model studied, that 
the diffusion coefficient, related to the position of the 
heavy impurity, is useless as a probe of the microscopic 
dynamics of the system to which it is coupled. On the 
contrary, we have presented evidence that the relaxation 
time of the MACF, related to the momentum of the heavy 
oscillator, shows a different behavior in going from the 
regular to the chaotic dynamical regime of the FPU chain 
as e increases. To the best of our knowledge this is the 
first time that such a manifestation of the microscopic 
dynamics of the FPU model can be observed in the be- 
havior of a physical macroscopic parameter characteristic 
of Brownian motion. So far we cannot generalize this re- 
sult to other models. This task will be left for future 
research. 
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